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ABSTRACT 

A multiaxial, isothermal, continuum damage mechanics model for fatigue of a 
unidirectional metal matrix composite volume element is presented. The model is 
phenomenological, stress based, and assumes a single scalar internal damage variable, the 
evolution of which is anisotropic. The development of the fatigue damage model, 
evolutionary law) is based on the definition ot an initially transversely isotropic fatigue 
limit surface, a static fracture surface, and a normalized stress amplitude function. The 
anisotropy of these surfaces and function, and therefore the model, is defined through 
physically meaningful invariants reflecting the local stress and material orientation. I his 
transversely isotropic model is shown, when taken to it’s isotropic limit, to directly 
simplify to a previously developed and validated isotropic fatigue continuum damage 

model. . . , 

Results of a nondimensional parametric study illustrate (1) the flexibility ot 
present formulation in attempting to characterize a class of composite materials and (2j 
the capability of the formulation in predicting anticipated qualitative trends in the 
fatigue behavior of unidirectional metal matrix composites. Also, specific material 
parameters representing an initial characterization of the composite system SiC/Ti 15-3 
and the matrix material (Ti 15-3) are reported. 


NOMENCLATURE 


Stresses 


"U 

S U 

^max 

Cauchy stress tensor 

deviatoric stress tensor 

normalized uniaxial maximum stress 

5 

normalized uniaxial mean stress 

a max 

uniaxial maximum stress 

a 

uniaxial mean stress 
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Invariants 


J 2 second invariant of the deviatoric stress (expressed for shear) 

Ij invariant representing transverse shear stress 
Ig invariant representing longitudinal shear stress 

Ig invariant representing maximum normal stress in the fiber direction 


Material Parameters 


cr u is the ultimate normal stress or static fracture stress 

<7^ normal stress endurance limit 

9t ratio of endurance limit to static fracture stress; 

M normalizing stress amplitude constant 

JC ratio of stress amplitude constant to static fracture stress; M/cr u 

< j 1 ratio of longitudinal to transverse ultimate normal stress 

ratio of longitudinal to transverse normal fatigue or endurance limit stress 
u/ m ratio of longitudinal to transverse normalizing normal stress amplitude 

77 u ratio of longitudinal to transverse shear static fracture stress 

rjy ratio of longitudinal to transverse shear fatigue limit stress 

7 ? m ratio of longitudinal to transverse shear normalizing stress amplitude 

a scaling factor for stress dependency in the fatigue damage exponent 
b,b’ parameters indicating effect of mean stress 
a fatigue damage variable exponent (can be a function of stress) 
exponent on the normalized stress amplitude 
6 angle between fibers and coordinate axis 

Subscripts 

( )p longitudinal properties 
( )rp transverse properties 
0 initial value 

Miscellaneous 

d- unit vector denoting local fiber direction 
D.j second order direction tensor; d-dj 
£.j Kronecker delta function 
Dp fatigue damage variable (measure) 

F^ fatigue limit surface 

F u static fracture surface 

F^ normalized stress amplitude function 
P combined hydrostatic mean stress measure 

£ parameter in P weighting the influence of the maximum hydrostatic stress 
N number of cycles 
Np number of cycles to failure 

< > Macauley brackets 
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INTRODUCTION 


In recent years, because of advances in scanning electron microscopy, acoustic 
emission techniques, and other nondestructive testing methods, a consensus has emerged 
that the nonlinear response of solids and ultimately their mechanical strength are 
dependent not only on the basic structure of the material but also on the type, 
distribution, size, and orientation of the defects in its structure. As a result, a 
relatively new branch of continuum mechanics, known as Continuum Damage Mechanics 
(CDM), has emerged. Continuum mechanics allows one to describe (over an appropriate 
representative volume element) the heterogeneous microprocesses involved during the 
straining of materials and structures at the macroscale. Similarly, CDM allows one to 
describe the material’s progressive deterioration (damage) from the virgin state (no 
damage) to the final state, corresponding generally to macrocrack initiation and 
propagation (or the ’’breaking up” of the representative volume element). Consequently, 
the life-limiting macroscopic properties such as rupture strength, fatigue life, or creep 
rupture lifetime of the material can be predicted. _ 

The tracking, or description, of the evolution of damage is accomplished by the 
introducing special thermodynamic (internal) field variables representing, in an 
appropriate statistical sense, the local distribution and density of defects. Numerous 
damage theories, both micromechanical [1-6] and phenomenological [7-11], have been 
proposed and discussed in the literature. There is great diversity in the mathematical 
nature of the damage variable(s) (e.g. , scalar [7,12,13], vectors [14,15] and tensors 
[12,16,17]), and thus the damage theories; this stems from the difficulty in directly 
measuring ’’damage” macroscopically and the degree of approximation with which the 
internal variables describe the salient aspects of the macroscopic effects of the 
microdefect kinematics. For a number of excellent review articles and books on this 
subject see references [7,12,18-24]. .... 

In this paper a phenomenological, isothermal, transversely isotropic, fatigue CDM 
model, with a scalar damage internal variable, applicable to unidirectional metallic 
composites, will be described. The present model is an extension of a previously 
validated CDM model developed at ONERA (Office Nationale d’Etudes et de Recherches 
Aerospatiales) [10,11,18,19,25,26,27,28] for isotropic monolithic metals. A brief review of 
the CDM work of Chaboche and his colleagues at ONERA is given in reference [291; it 
provides a foundation for the present extension of the ONERA models to transversely 
isotropic materials (e.g., hexagonally packed metal matrix composites). 


A CDM FATIGUE MODEL FOR INITIALLY ANISOTROPIC MATERIALS 

Different damage variables are associated with different damage processes, such as 
creep, fatigue, ductile, and brittle damage, and have been discussed^ in the literature 
extensively. Here we will confine ourselves to the modeling of fatigue damage by using 
a scalar damage measure. Figure 1 clearly illustrates the particular damage mechanisms 
considered in this study of the fatigue of monolithic metals. In Fig. 1, the fatigue crack 
initiation and growth process is illustrated schematically; the damage (Dp) is associated, 

macroscopically, with the initiation and propagation of transgranular defects (e g- , slip 
bands and microcracks). The possible equivalence between the definition of Dp by the 

effective stress concept, the definition of Dp in terms of the remaining life concept, and 

the quantification of physical damage in terms of micro-cracking has been shown by 
microcrack measurements made by Cailletaud and Levaillant [30], Hua and Socie [31], 
and Socie et al [32] 

The Continuum Damage Mechanics approach, which supports the previously 
described global measure of fatigue damage in monolithic isotropic materials, presents 
two theoretical deficiencies which make it inconsistent with a continuum approach: (1) 

the surface character of fatigue damage and (2) the small number of large defects 
present during the propagation stage. Interestingly enough however, these deficiencies do 
not appear to affect its ability to correctly describe many experimental results 
(performed under uniaxial loading) in terms of cumulative fatigue damage. 
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Figure t. —Different damage modes and the associated scale in a monolithic metal. 


We believe that these two theoretical deficiencies can be removed by establishing a 
representative volume element (RVE) on the mesostructural scale of a composite 
material; then the surface character of fatigue will transfer to the interface region of 
each constituent (or phase), and the number of large defects occurring during the 
propagation period will similarly increase. This hypothesis was borne out by observation 
of experiments conducted under specific thermal and mechanical loading histories [33,34]. 
There are, however, a number of surface and environmental effects which will not 
necessarily transfer to the fibers, for example those effects that are associated with 
oxidation in the SiC/Ti 15-3 system [35]. A nonrigorous definition of an RVE for a 
given composite material is given in Fig. 2.. In short, it requires that the RVE include 
a sufficient number of unit cells to allow (statistically) the heterogeneous nature of the 
material to be homogenized. 

Figure 3 is the extension of Fig. 1 to the mesoscale; it illustrates the 
mesostructural damage (i.e. , the microcracks and mesocracks) that might result when a 
composite material is subjected to fatigue loading. Here, damage mechanisms similar to 
those on the microstructural scale are postulated to occur on the mesostructural scale; 
the surface or interface of a constituent (the fiber) plays the role (on the mesostructural 
scale! of the grain boundary of a constituent (grain) on the microscopic scale. 
Furthermore, because of the internal structure of the material, a sufficient number of 
defects will be present, thus allowing a theoretically consistent continuum representation 
of fatigue damage to be constructed. 

In addition to matrix cracking and the formation of microcracks around the fibers, 
other damage modes must be included, such as fiber breaking, fiber/matrix debonding, 
and interlaminar cracking. It has been noted that damage occurs in the form of 
different multiple cracking modes and that there is no isolated single crack that 
dominates the development of damage. Similarly in composite laminates, damage develops 
along preferred orientations; for example, matrix cracks in off-axis plies typically are 
channeled by the fibers in those plies and interlaminar planar cracks grow along fibers in 
the neighboring plies [36]. 
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Figure 2.— Nonrigorous definition of representative vol- 
ume element (RVE); s u * displacement field, s a = 
stress field. 
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Fatigue damage: initiation and propagation of 
transgranular defects. 

Figure 3.— Different damage modes and the associated scale in a metallic composite. 
Fibers on mesostructural scale are analogous to grain boundaries on microstructural 
scale. 


These direction-oriented damage modes suggest the need for a direction-dependent 
damage variable (e.g., a vector or tensor). Here, however, on the basis of the strong 
initial anisotropy of the composite, we assumed that the damage measure is a scalar, 
with the evolution of damage being anisotropic and associated with the preferred direction 
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of the material. For example, if the composite is loaded in the fiber direction, the 
evolution of damage would be small relative to the case in which the composite is 
loaded normal to the fiber. However, because of the assumed scalar measure, the 
distribution of damage would be identical in all directions. Thus the accumulation but 
not the magnitude of damage is direction dependent. The approach taken to model this 
damage was the stress based approach detailed in references [12] and 129], wherein the 
definition of effective stress and the concept of remaining life were utilized to measure 
and interpret damage. 

Multiaxial Statement 

The ONERA model was extended by introducing anisotropic damage surfaces with 
appropriate stress invariants [37-39] representing stress states that are likely to strongly 
influence the various damage modes in metallic composites. For instance, we assumed 
(1) that the transverse shear stress (1^) and thus implicitly the transverse normal stress, 

dictates matrix cracking; (2) that the longitudinal shear stress (Ig) dictates interfacial 

degradation; and (3) that the maximum normal stress (1^) in the fiber direction dictates 

fiber breakage. In this way, we believe that the following fatigue limit and static 
fracture surfaces and the normalized stress amplitude function contain the most pertinent 
damage-producing stress measures. Let 


F 0 - 


1 

T77 


(4 & f)~i)i 1 + Tf 


( 4 W 1 ) 


() 


^ + 


where the fatigue limit surface 1 is defined as 


(i) 


i 


F « - 1 = 2 m f c F (fO lcr ij (t) " a i “ 1 


( 2 ) 


the static fracture surface is 


u ‘ T“ ‘ 


1 - F„ = 1 - max F/ 


and the normalized stress amplitude is; 


( 3 ) 


1 


F m = F m( ff ij) = 2 m f m f " 'ijW 


( 4 ) 


The subscript ( ) in Eq.(l) is to be replaced by either fl, u, or m. Material parameters 
(some of them varying with the mean stress) that describe the initial anisotropy of the 
surfaces are defined as follows: 


U L ~ a {l L ’ 

= T uJ T 1Lj 

w il - 


m L = M L ; 

= M L/ M T : 


U L = \ 

In = T »J T , 


V U T 


w u = u L /u T 


*The authors acknowledge that not all monolithic materials exhibit a mar ked fatigue or 
endurance limit fe.g., nonferrous alloys); however we assume (on the basis of 
experimental evidence [40]) that a sufficiently large class of composite materials will 
possess such a limit. Furthermore, we assert that fatigue lives beyond some specified 

n 

value, say 10 , can be considered infinite for most practical design purposes and thereby 
allow the assumption of an endurance limit. 
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Note that the subscripts L and T denote, respectively, the longitudinal (along the fiber) 
and transverse (normal to the fiber) directions. 

The physically meaningful invariants Ip I 2 , and Ig, respectively, are defined as (cf. 


[39,41]); 


J 1 ~ J 2 ” 1 + 4 T 3 

i 2 = i - i 3 


( 5 ) 

( 6 ) 
( 7 ) 


in which 


D; ; 


S U 


“ 2 S ij S ij 
= D ij 

= D ij S jk S ki 
i 

= a ij ~ T °kk 5 ij 


( 8 ) 


and d- (i = 1, 2, 3) represents the components of a unit vector denoting the local fiber 
direction. 

Now, from the fatigue limit, and static fracture surfaces and normalized stress 
amplitude function we can formulate the muitiaxial representation of the extended 
fatigue damage evolutionary law; this is valid for materials that possess initial 
transversely isotropic material symmetry (e.g., hexagonally packed unidirectional 
composite materials): 


dD F = 


1 - (1-Dp) 


£+1 


1 a 

r 


m 

. 

( 1 Dp ) 


dN 


( 9 ) 




A 

7 

w- 

V 


with 

a = 1 

- a - — - — — 

(10) 



< 1 -F u > 


and 

U L 

<J {L (0) (1 - 3b’ P) 
U L 

(11) 


II 

S 

M q (1 - 3b P) 

(12) 


where a, b, b\ (3 , (0), a , M n are material constants and henceforth, 

u h U L U L 


< X > 


= 0, X < 0 

[ = x, X > 0 


defines the Macauley bracket. 

One additional modification, besides the inclusion of a preferred direction, has been 
introduced into the preceding fatigue damage model: a hydrostatic stress measure P that 
combines both the Sines [42] and Crossland [43] criteria, that is, 

P = C a Hmax + f 1 “ 0 °Havg. ^ 13) 
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where 


A 

ff Havg = 2 { m * x [ <T ii( t )] +m j n Ki^W 
ff Hmax = 


(14) 

(15) 


are the mean hydrostatic stress and maximum hydrostatic stress, respectively. 

We feel this alternate hydrostatic measure (P) is required, since a judgment about 
the most applicable criterion is unavailable because of a lack of experimental data on 
metallic composites. In the ONERA model, only the Sines 1 criterion was considered 
(i.e., ( = 0). Finally, if we assume that the ratios <j u , t? u , w m , and rj m , are 

each equal to one then the isotropic ONERA model is easily recovered with only a 
slight modification to the coefficient a in equation (10). 

Nondimensional Uniaxial Simplification 

In order to conveniently examine the behavior of a variety of composite systems, 
consider a uniaxial stress state a, fibers oriented at an angle 0 from the x-axis (i.e., 
dj = (cos 6 t sin 0 , 0) in the x-y plane), and the introduction of the following 

dimensionless parameters. Here both the maximum and mean stresses have been 
normalized relative to the ultimate (or static fracture) stress (cr^ ) in the fiber direction; 

L 

that is, 


and 


max 


max 


(16) 


5 = 


respectively. Thus, the uniaxial simplification of equations (9) and (10) are 


dDj 

dN 




and 


f l ( S max~ ? )' 1 ^ 


a = 1 - a 


i — yw” s ] 

~ u max J 


where 


(4^-1) 


A 1 + — J A 2 + A 3> 

1 ( 4 ‘ J ff 1 ) 

= & A i + A 2 + V 

( 4<-l) 

= {(4<-l) A 1 + A 2 + A 3 ) 


(17) 


(18) 


(19) 
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M = 


M, 


9t = 


u L 
'“l 

r 

Hi 


( 20 ) 


( 21 ) 


A, = j — g cos 2 , — g sin 2 , 4- gg (2 cos , — sin ,) 

A 2 = g cos 2 , 4- g sin 2 , - g (2 cos 2 , - sin 2 ,) 2 
A 3 = \ (2 cos 2 , - sin 2 ,) 2 

Note that M L (P), u^P), <r u (P), and the _w and i) ratios are independent of cr H since 

L 

the modification factor cancels. Again, if u/q and tjq are equal to 1 (i.e., the material 
is isotropic) then equations (17) and (18) become 


dD p = [l-(l-D F )^ +1 ]“ 


S -7 
max 


and 

with 

and 


a = 1 - a 


M l (1-Dp) 

. ( S max -S f<^) 

(^u~^max) 

S f /5) = S f /0)+[l - b S f /0)] 5 

a ‘ a fe) 


0 


dN 


( 22 ) 


(23) 


which is identical in form to the uniaxial ONERA model I (see [29]). 

Integrating equation (17) for N cycles, where S max and R are held fixed, leads to 


N = 


[l-(l-D F )' J+1 ] (1 a) r 1-0 

(T+?) ( 1— a”) [y ® (S> max b)i 


(24) 


By cycling the material to failure (Np), which is defined to occur when Dp — 1, we 
find that 


N F - (i+/»)(i-a) 


JWW (S - 3) 

v m v max } \ 


-p 


(25) 
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or 


[ tJ ( S m a x *5") ] ^ I [ 1-y/ # u S^ax] \ ^26) 

^ ” UTP) a ( Sa,ax-5)-lJ / 

An expression describing the accumulation of damage (for a given loading cycle) in 
terms of the remaining life ratio can be obtained by substituting equation (25) into 
equation (24) and rewriting the resulting equation: 

D p = 1 - [l - (N/Np) 1 ^ 1- ^] 1 ^ 1 ^ (27) 

Note that the accumulation of damage is nonlinear and is also a function of fiber 
orientation , as is the number of cycles to failure. 

Similarly, if a two-level loading test is considered, one obtains, by integrating 
equation (17) in two steps, the same expression as in the case of isotropic materials: 

N 2 /N f = 1 - (N^Np )P p = (1 - a 2 )/( 1 - Oj) (28) 

2 X 

Here, though a., a 9 , N p , and N P all depend on the orientation of the fibers (or 
1 1 *1 f 2 

conversely the load). Note that N P and N P are the failure lives (on the S~N curve) 

• 1 *2 

for the first and second loading conditions, respectively; whereas N^ is the number of 
cycles applied at the first loading level and Ng is the remaining life at the second level. 


RESULTS: PARAMETRIC STUDY 


Direction Independent Parameters 


For an isotropic material (i.e., Wq = tj,\ = 1) the effect 


on the S-N curve of 


varying the parameters a, and the mean stress (5) was examined, see reference 

[29J. To simplify and yet not limit the study, the parameters <7^ and were taken 

to be independent of hydrostatic stress; that is to say, b and b' (given in equations (11) 
and (12)) are zero. The isotropic baseline values ^ assumed were those corresponding to 
304 stainless steel [12] (see Table 1). Because material isotropy was assumed, equation 
(26) becomes 


N F - (l+0)a 


a x' 

—fT 


/ [ 1 - Smax ] 

** '[S m ax-(3+ &) y 


(29) 


Equation (29) suggests that if either the mean stress (5) or fatigue ratio ( #) are 
modified, the load level at which an endurance limit is reached is also changed; that is, 
as either value is increased so is the endurance level. Furthermore, equation (29) 
indicates that translation of the S-N curve can be achieved by varying the parameters 
JC and a. Increasing JC has been shown (see [29]) to shift the S-N curve to the right 
(increase the number of cycles to failure at a given load), whereas increasing a shifted 
the S-N curve to the left (decreased the number of cycles to failure at a given load). 

Thus only the product (a JC~^), not the individual values of a and JC^ is important in 
determining the life at a given load level. However, the actual value of a does 
strongly affect the damage accumulation, and therefore iL_a ny other type of damage ($•£• 
creep) is present , both parameters JC and a need to be determined explicitly. 
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Table 1 : Dimensionless Parameters for Isotropic and Anisotropic Materials 


Parameters 

Values for Baseline 
Transversely Isotropic 

Values for 
SiC/Ti 15-3 

Values for 
Ti 15-3 

5e 

0.25 

0.18 

0.16 

JC 

4.5 

14.42 

7.03 

P 

3.0 

1.842 

2.27 

a 

0.1 

0.065 

0.230 

5 

0.0 

S/2 

S/2 


5.0 

11.8 

1.0 

U u 

5.0 

12.482 

1.0 

w n 

5.0 

5.5 

1.0 

V 

2.0 

1.0 

1.0 


2.0 

1.0 

1.0 

% 

2.0 

1.0 

1.0 

\ 


225. 

128. 


Finally, from equation (29) one can deduce that perturbations in the parameter 0 will 
both translate the S-N curve and affect the abruptness of the transition from infinite life 
to immediate failure. Clearly, 0 = 0 is the lower limit and will give the most abrupt 
S-N curve for a fixed set of material parameters [29]. 

Parameters Defining Initial Transverse Isotropy 

Assuming the baseline transversely isotropic material parameters given in Table 1, 
we examine the modification of the S-N curve with fiber orientation (or equivalently, 
load orientation relative to a fixed fiber orientation, see Fig. 4). As might be expected 
(since a composite material is designed to be stronger when loaded parallel to the fiber 
direction), both the static fracture strength and fatigue limit decreased with an increase 
in the Coer orientation angle. The reduction in load-carrying ability is clearly 
dependent on the 11 degree" of anisotropy and is represented in this fatigue damage model 
by the u) and rj ratios in the ultimate stress ( fatigue limit ( and normalizing 


Fiber 



Figure 4.— Illustration of the degradation in load-life response 
with variation in fiber orientation 0. 
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stress amplitude ( && m ) (see Eq. (26)). Here, with the given baseline material 

parameter set, all three anisotropic functions are considered to have the same degree of 
initial anisotropy. In this section an examination of the impact of perturbing the 
various anisotropic measures will be undertaken. 

Equations (19) to ( 21 ) and (26) indicate that the angle dependency is manifested 
through the functions A p A 2 , ana A^, which are associated with the longitudinal shear, 

transverse shear, and normal stress components, respectively. This angle dependency is 
clearly shown in Fig. 5 . Furthermore, Fig. 5 indicates that when 0= 0 # , any changes in 
the measures of anisotropy (i.e., u or rj) will have no effect; this suggests that initial 
characterization of material parameters a, b, /3, and JC, as discussed in [29], be 
conducted with longitudinally reinforced and loaded specimens. Similarly, when 0 = 90% 
any change in the shear "strength” measures (i.e., r? u , 77 ^, and rf m ) will not have any 

impact in the resulting S-N curves; this suggests that independent evaluation of the oj , 

Wy, and can be achieved by comparing the longitudinal and transverse specimen 

response. 



Figure 5. — The angle dependency of the longitudinal shear 
(A-j), transverse shear (A 2 ) and normal stress (A 3 ) compo- 
nents. 


Variation of u . 

-u J 


(Jr/ and 
-it m 


Figure 6 illustrates the effect on the S-N curve (corresponding to a fiber 
orientation of 8 = 15°) of varying the ultimate strength ratio (u/ u =< 7 u f a n )• The 

L T 


values of a-% are taken to be 2 , 4 , 8 , and 16, while all other parameters are held fixed 

at the baseline values given in Table 1 . As might be expected increasing the ultimate 
strength ratio affects only the low cycle fatigue portion of the curve (by decreasing the 
static fracture stress); the endurance limit remains unaffected. Conversely, increasing the 
fatigue limit ratio (i.e., = 2 , 4 , 8 , and 16) and holding all other parameters fixed, as 

shown in Fig. 7, decreases the endurance limit and leaves the ultimate strength 
unchanged. Note, that although the present model provides significant flexibility, in that 
uj u and can be varied independently, in reality a relationship between the ultimate 

strength and the endurance limit would be expected, thus placing restrictions on the 
ranges of these parameters. 
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By varying the normalizing stress ratio u/ m , a horizontal translation of the 

transition portion of the S-N curve can be obtained while the ultimate stress and fatigue 
limit remain unchanged. An example is shown in Fig. 8, where takes on the values 

of 2, 4, 8, and 16, the fiber orientation is 15* , and all other parameters again remain 
fixed. Clearly, increasing the ratio shifts the S-N curve to the left thereby decreasing 
the number of cycles to failure for a given load. Thus, by merely assuming different 
degrees of anisotropy for the three functions and (1) either end of the 

S-N curve can be modified while the other end remains unchanged, or (2) the center 
portion can be horizontally translated, while the ultimate and fatigue limit stresses 
remain unchanged. This suggests that the model has sufficient flexibility to fit a broad 
class of materials. 



Figure 8.—Effect on the S-N curve of varying the ratio a> m at 
an angle of 15° while holding all other parameters fixed at 
the baseline values. 


The impact on the S-N curve of varying the degree of anisotropy, for example u/ u , 

can be seen for all fiber orientation angles by plotting the functions , 

and yj # In Fig. 9 only is shown versus angle of orientation 0 } since the 

other two functions behave similarly, except for being scaled by a multiplying factor. It 

is evident in Fig. 9 that for a specified ratio a; the function increases as the 

angle 0 increases; however, there is an angle, depending upon the magnitude of w , at 

which this rate of increase decreases. This explains the incrementally smaller changes in 
the response curve with increasingly larger angles (i.e., variation in the S-N curve, see 
Fig. 4). Finally, another important point that can be discerned from Fig. 9, is the need 
for u to always be greater than or equal to r?, since if u < r\ y a minimum occurs in the 
function at an angle other than zero, which violates physical reasoning. 
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Figure 9. — Impact of varying w u when t] u = 2.0 and all other 
parameters are those in Table I. 


Variation of t? u , ij f f , and rj m 

Variation of the shear ratios 77 , 7/^ and r) m has an impact on the overall trends, 

but not actual magnitudes, of the S-N curve, similar to variations in the respective 
normal stress ratios. Examples involving the variations of each ratio are given in [29]. 
The interaction between the shear strength and the normal strength ratio is shown in 
Fig. 10, for 5.0, and tj u = 1, 1.5, 2, 3, 4, and 5; and in Fig. 11, for w u =16.0 and 

77 u = 1.0, 1.5, 2, 4, 8, and 16. Note how in both figures the shape of the function 

is greatly affected, particularly for low values of 77. Also it is clear that when 77 = 

1, a maximum at some intermediate angle is obtained. The ramification of this 
intermediate maximum is unknown at this time. As with the cases involving the 
normal stress ratios and Wy , one would expect a relationship to exist between rj n 

and rjy such that the practical ranges of these parameters are not independent. 



Figure 1 0. — Impact of varying ti u when <u u = 5.0 and all other 
parameters are those in Table I. 
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Figure 11 . — Impact of varying ti u when oj u = 16.0 and all other 
parameters are those in Table I. 


CHARACTERIZATION EXPERIMENTS 

The required exploration and characterization experiments for the present model are 
completely described in [29] . Here, we want to briefly touch on the testing required, 
both in tension and torsion in the longitudinal and transverse directions, to define the 
extent of initial anisotropy (i.e., the ratios a> u , t? u> rj^ a> m , and rj m ) of a 

hexagonally packed unidirectional composite material. Figure 12 illustrates schematically 
the expected S-N curves for a composite material loaded longitudinally (9 = 0 8 ) and 
transversely (9 = 90 8 ) at a given mean stress. From such data, as well as the 
associated shear S-N curves, the various uj and tj ratios can be determined. Note that 
in the present theory, 0, a, and b are assumed to be independent of fiber orientation 
(i.e., d*;. This assumption may need to be modified if experimental evidence suggests 

otherwise. 

Clearly, the four S-N curves (Fig. 12, repeated once for tension and once for 
shear) as well as other experiments discussed in [29] suggest a significant, and poten 
tially quite expensive, experimental program. This experimental program may, however, 
be able to be reduced and augmented through numerical simulation by using a suitable 
homogenization technique with periodic boundary conditions. Furthermore, utilization oi 
a homogenization technique should allow micromechanical effects (e.g., bond strength, 
volume fraction, etc.) to be included in this phenomenological theory, through the 
parameters a , o> u , ?? u , uty and 

Recently, a model metal matrix composite system (i.e., SiC/Ti 15-3) and. matrix 
material (i.e., Ti 15—3) were used in the correlation of the present fatigue CDM model. 
Figure 13 shows both the correlations (lines) and the experimental data (symbols) for a 
35 fiber %, longitudinal and transverse composite and the Ti 15—3 matrix material at 
427° C. Table 1 lists the specific individual material parameters used. As suggested 

earlier, the longitudinal response was used to obtain the constants a, /?, a , and 

Ju 

a r , . Then while these constants were held fixed, the three transverse data points 

were fit to obtain the "strengths" of anisotropy, and Note that the shear 

ratios, 77 , 7 7^ and r\ m were all assumed to be 1.0 (i.e., isotropic in shear) and the 
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effect of mean stress on the fatigue limit and stress amplitude (M) was neglected (i.e., 
b=Q) because of the lack of data. Finally, in Fig. 13, the solid square symbol 
associated with the 90* response was not used in the correlation data set for the 
composite response and therefore represents a check on the predictive capability of the 
present model at the coupon level. The accuracy is within a factor of 2. 



0 1 2 3 4 5 


log (N) 

Figure 1 2.— Expected S-N curves for a unidirectional 
metallic composite loaded longitudinally (9 = 0) and 
transversely (6 = 90). 
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Figure 13. — Correlation of CDM model with SiC/Ti-15-3 com- 
posite and Ti-1 5-3 matrix coupon data. 
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CONCLUSIONS 


An isothermal, multiaxial, continuum damage mechanics model for fatigue of a 
unidirectional metal matrix composite volume element has been presented. The intended 
applications are reinforced structures in which the fiber direction may vary throughout 
but in which a single fiber direction can be identified locally (local transverse isotropy) 
within a given volume element. The model is phenomenological, stress based, and 
assumes a single scalar internal damage variable, whose evolution is anisotropic and is 
associated with the initiation and propagation of transgranular defects. 

The development of this model (i.e. , the fatigue damage evolutionary law) is 
founded on the definition of an initially transversely isotropic fatigue limit surface, static 
fracture surface and a normalized stress amplitude function. The anisotropy of the 
model is defined through physically meaningful invariants reflecting the local stress and 
material orientation. The model has been shown, when taken to its isotropic limit, to 
directly simplify to a previously developed and validated fatigue continuum damage 
model. 

Results of a nondimensional parametric study illustrate (1) the flexibility of the 
present formulation in attempts to characterize a large class of composite materials and 
(2) the ability of the formulation to predict anticipated qualitative trends in the fatigue 
behavior of unidirectional metal matrix composites. Additionally, an example correlation 
of a metal matrix composite system (i.e., SiC/Ti 15—3) and metal matrix material (i.e., 
Ti 15-3) was shown. 

Two potential drawbacks to the present formulation are (1) the scalar damage 
measure employed and (2) the potentially expensive experimental program required. 
Future work in this area will include an examination of the present formulation’s 
usefulness in predicting the fatigue life of high-temperature applications; initiation of an 
experimental characterization and validation program; and the numerical implementation 
into a post processing life prediction computer code. 
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